The magnetic form factor of the deuteron in chiral effective field theory 
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We calculate the magnetic form factor of the deuteron up to 0{eP*) in the chiral EFT expansion of 
the electromagnetic current operator. The two LECs which enter the two-body part of the isoscalar 
NN three-current operator are fit to experimental data, and the resulting values are of natural size. 
The 0{eP'^) description of Gm agrees with data for momentum transfers < 0.35 GeV^. 
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Introduction: In the past two decades chiral effective 
field theory (xEFT) was fruitfully applied to few-nucleon 
dynamics (see Refs. [H, 0| for recent reviews). Two- 
nucleon potentials at next-to-next-to-next-to-leading or- 
der (N^LO) in the chiral expansion were developed [1, Q 
which accurately describe low-energy scattering data and 
the static properties of the deuteron. Higher-order cor- 
rections to the three-nucleon force are presently under 
investigation, see e.g. [l^l, although discussions regard- 
ing non-perturbative renormalization of the Schrodinger 
equation and implications for the xEFT power counting 
continue, see [5|-l9| for samples of different views. In par- 
allel to these developments in the strong sector, much ef- 
fort has been devoted to pionic and electroweak reactions 
in few-nucleon systems, see [TT| - [l3j for recent examples. 



Electromagnetic reactions on light nuclei such as elas- 
tic electron scattering, photo-/electrodisintegration and 
radiative capture have been extensively studied in nu- 
clear physics. In the single-photon approximation, their 
theoretical description requires knowledge of the electro- 
magnetic current operator, which should be constructed 
consistently with the nuclear Hamiltonian. The deriva- 
tion of exchange currents in xEFT was first addressed in 
the seminal paper by Park et al., [13], who, however, lim- 
ited themselves to threshold kinematics |q| ^ M.^ with 
denoting the pion mass. Recently, this work was ex- 
tended to the general kinematics suitable to study, e.g., 
electron scattering off light nuclei at mometum transfer 
of |q| of order M^r by the JLab-Pisa [T5l - [l7| and Bochum- 
Bonn groups [H, [3- Here and in what follows, we dis- 
cuss the expansion of the irreducible two-nucleon opera- 
tors Jo and J in powers of P = {p, TOtt ) / A with A denoting 
the hard scale in the theory, e.g. the cutoff (~ 600 MeV) 
used in calculations. In this expansion the leading-loop 
order is eP^. However, most of the corrections to the two- 
body pieces of the two-nucleon current and charge oper- 
ators at this order are of isovector type and thus do not 
contribute to the deuteron form factors. In particular, 
up to this order, the only two-body contributions to the 
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isoscalar charge density operator, ' emerge from the 
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leading relativistic corrections of one-pion range so that 

(s) 

Jg is parameter free. The impact of these corrections 
on the deuteron charge and qua drup ole form factors, Gc 
and Gq is studied in Refs. [20, UM- I'^ these works the 
deuteron wave functions obtained from xEFT potentials 
at various orders were used to compute Gc and Gq (see 
also Refs. [l^jH^I for earlier work along the same lines). 
Good agreement with the compilation of elastic electron- 
deuteron data from Ref. [22| was then found for both 
form factors in the kinematic range < 0.35 GeV^, 
provided factorization was employed in order to account 
for single-nucleon structure. 

On the other hand, the isoscalar two-nucleon current 
operator, J*^*' has two two-body contributions at order 
eP^: one from a short- distance operator and one of one- 
pion range. The impact of these terms on the magnetic 
moments of the deuteron and trinucleons was examined 
in Ref. [251] . However, markedly more information on the 
interplay of these terms with each other and with one- 
body mechanisms is available via the |qp-dependence of 
observables. In this work we present a study in this 
direction, using xEFT expressions for J*-") derived in 
18 . 19 1 to extend the predictions given for Gm in 
23| to 0{eP'^). We discuss the relevant terms 
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in the current operator and use the data on the magnetic 
form factor of the deuteron at low values of |qp to deter- 
mine two unknown low-energy constants (LECs). The 
0{eP'^) xEFT results thereby obtained accurately de- 
scribe experimental data on Gm in the kinematic range 
Q2 < 0.35 GcV^ This, together with the findings of 
Ref. 21], provides a full set of results for elastic electron- 
deuteron form factors at 0{eP'^). 

In the next section we describe the anatomy of the 
calculation and outline the relevant terms in the two- 
nucleon current operator. This is followed by a discussion 
of our results, including those for the LECs. We finish 
by summarizing. 

Anatomy of the calculation: The magnetic form factor 
of the deuteron we are focused in this work is related to 
the Breit-frame matrix element of the four-current oper- 
ator Jfj_ according to the well-known relation 
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where = + i,P and 77 = ]q]^/(4TO^) with q = 
Pg — Pe denoting the momentum transfer and rrid the 
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deuteron mass. (Since we work in the Breit frame we 
have = |qp.) The deuteron states are labelled by 
the projection of its spin along the direction of q. Both 
the deuteron wave functions and the current operator 
appearing on the right-hand side of the above equation 
are calculated order-by-order in xEFT. 

We now briefly describe the xEFT expansion of the 
two-nucleon current operator, J^, as it pertains to the 
calculation of the deuteron form factors. We employ 
Weinberg's power counting throughout this work which 
makes use of naive dimensional analysis to determine the 
significance of various contributions. The leading con- 
tribution to the charge density, Jq, is given by an 
photon coupling to a point proton at order e. Nucleon- 
structure corrections start contributing to the one-body 
current at order eP^ [next-to-leading order (NLO)]. The 
first isoscalar two-body contribution is generated from 
a tree- level pion-exchange diagram at order eP^, pro- 
vided the nucleon mass is counted as mA? ~ A^/M^ ^ A 
[3) [3 ^ • This correction is associated with a relativis- 
tic correction to the one-pion-exchange part of the NN 
potential. There are numerous other corrections to the 
two-body part of the charge operator at leading-loop or- 
der, eP^, from one- and two-pion exchange diagrams and 
from pion loops involving the lowest-order contact inter- 
actions, but this is the only isoscalar effect. The explicit 
form of all terms can be found in Refs. [Tsl. ITgj. 

The chiral expansion of the three-current starts at or- 
der eP with the single-nucleon contributions. The first 
two-body terms emerge from tree-level one-pion exchange 
diagrams at order eP^ [NLO]. The next two-body correc- 
tions to J occur at order eP* from pion loops and tree di- 
agrams involving higher-order vertices from the effective 
Lagrangian. The two-pion exchange contributions are 
parameter- free [l8| . while the one-pion exchange terms 
depend on the LECs d^, dg, dig, d2i and (^22 entering 
jC^^ [13, HHHl ■ However, the only long-range two-body 
mechanism in J'^'*-' at this order is proportional to dg. 
While this LEG could, in principle, be constrained by 
pion photoproduction data, in p ractice these data pro- 
vide little information on dg j29l . Isoj . 

Gm is of particular interest at 0{eP'^) because it is 
there that the first short-distance NN physics not de- 
termined by NN scattering and gauge invariance ap- 
pears. This is represented by the simplest Ml isoscalar 
four-nucleon-one-photon contact term in the xEFT La- 
grangian, which is of the form [19., 23^, ^] : 
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The low-energy constant (LEG) L2 that appears in 
Eq. ([2]) must be extracted from data on electromagnetic 
reactions in the two-nucleon system. 



^ In the nomenclature of Refs. Il9l|. these contributions appear 
at O(eQ). 



The combination of these two effects yields a two-body 
isoscalar current operator J^^^ (Tsl. [Toj: 



J(s) o " 1 " ^ r 1 

IB = 2e^d9 • T2 ^2^:^ [qi X q] 
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ieL2 (0-1 + 0-2) X qi + (1 2) , (3) 

where q labels the photon momentum and qi/2 labels 
the momentum transfer on nucleon one/two respectively. 
Since G m is determined completely by the one-body part 
of J^*-* up to 0(eP'^) the total form factor is thus 
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Here we use factorization to compute Jib, i.e. we write: 
^[f = i^[4^)(g2)2p+ + zGif (g2)(ai X q)+], (5) 

with p the momentum of the struck nucleon, and 

and G\j the isoscalar single-nucleon form factors, for 
which we take the parameterization of Ref. [s^ . The use 

of this ansatz for the one-body part of J^*-* is equivalent 
(up to corrections that begin only two orders beyond the 
order to which we work) to making a yEFT expansion for 
the "body" form factors Dm and De This allows us 
to focus on the momentum transfer at which the xEFT 
expansion for the NN current operator J breaks down, 
without having to worry whether the theory is doing a 
good job of describing isoscalar nucleon structure. 

Results: We now evaluate the matrix elements in 
Eq. ([U with a variety of xEFT deuteron wave functions 
computed with the NLO and NNLO xEFT potentials 
and different values of the cutoffs A in the Lippmann- 
Schwinger equation and A in the spectral function. The 
result found for Gm with LO xEFT wave functions and 
the leading piece of J^*\ denoted here as 0{eP), was 
computed in Ref. [s^l- Gorrections to this come both 
from higher-order pieces of the NN potential, V, which 
affect the wave function, and from the corrections to J*^''-' 
discussed in the previous section. The NNLO xEFT po- 
tential includes all effects up to O(P^) relative to leading 
(in this counting), so its deuteron wave function, when 
combined with the O(eP^) J^''^ yields a xEFT calcula- 
tion for Gm which includes all effects up to 0{eP'^). 

The pertinent matrix elements are computed via 
Monte-Garlo (MC) integration. To increase efficiency, 
we use importance sampling with the weight function of 
Ref. j35|]: 



p(k) =p(/c) 



(r-3)(r-2)(r-l) C 



r-3 



(fc + cy 



(6) 



The functional form oip{k) is chosen such that the weight 
function is maximal at the origin, reflecting the large 
S'-wave component of the deuteron wave function. The 
parameters C and r control the vanishing of the weight 
function at large momenta and are tuned to optimal val- 
ues (in terms of the efficiency of the MG integration) 
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FIG. 1: The magnetic formfactor Gm as a function of |q|. Ex- 
perimental data for the magnetic moment is from [SQ]. The 
remaining data are from the parameterization of 13711 (upward 
triangles) and scattering experiments reported in (down- 
ward triangles), 39] (squares) and (solid dots). The light 
blue (dark red) band represents the results with NLO (NNLO) 
wave functions, and j'*' computed up to 0(eP*). 



by calculating the expectation value of the one-pion ex- 
change potential yielding C = 1 GeV and r — 

As in Ref. [HI we perform a path average over sev- 
eral runs. We use 2730 sample points and and the path 
average is performed for 3000 runs. Analysis of the run- 
to-run fluctuations indicates a final answer with better 
than 1% precision throughout the momentum range of 
— 800 MeV. At several points we compared this MC 
answer to calculations using quadrature methods, and 
always found agreement within the precision claimed. 

We adopt the following procedure to determine the val- 
ues of the two LECs entering J^*'. First, we fix the value 
of L2 for a given dg by demanding that the magnetic mo- 
ment of the deuteron is reproduced. We then perform a 
X^-fit to the experimental data for |q| < 400 MeV (in- 
cluding four points from the parametrization of Ref. (37[ ) 
to determine c?g . Our attempts to use even lower- |q| data 
for this fit resulted in unstable answers, reflecting the in- 
sensitivity of Gm to this LEG at small values of |q|. 

The results of this procedure are shown in Fig. [T] 
The light/blue (dark/red) band is obtained using wave 
functions computed with the NLO (NNLO) xEFT po- 
tential. The width of the band shows the variation of 
the prediction as A and A are changed in the range 
A = 400 .. . 550 MeV (A = 450 .. . 600 MeV) at NLO 
(NNLO) and A = 500 .. . 700 MeV. The cutoff varia- 
tion is reduced at NNLO, and the data well described 
for g2 < 0.35 GeV2. 

In order to assess the momentum transfer at which 
the xEFT expansion for J*^"^ breaks down, in Fig. [5] we 
show the size of different contributions to the final result. 
This time the bands represent the impulse approxima- 



FIG. 2: The magnetic formfactor Gm as a function of |q|. 
The light blue (dark red) band represents the results with 
NLO (NNLO) wave functions using the impulse approxima- 
tion. The dashed/dotted lines are the contributions from two- 
body pieces of j'"', as described in the text. For remaining 
notation see Fig. [T] 

tion resuh obtained with NLO (light /blue) and NNLO 
(dark/red) wave functions. The dotted (dashed) line is 

(s) 

the effect from the piece of 3\g that is proportional to 
L2 (dg). For both two-body matrix elements, we show 
results averaged over the five cutoff combinations con- 
sidered, with the light blue lines showing the NLO case, 
and the dark red lines obtained with NNLO wave func- 
tions. We estimate the breakdown scale of the EFT ex- 
pansion by values of momentum transfer at which the 
0{eP'^) two-body contributions start becoming compa- 
rable to the effect of the 0{eP) (impulse-approximation) 
piece of the current. Fig. [2] shows that the smaller two- 
body contributions to Gm found with the NNLO wave 
function delay the breakdown of the expansion. Even 
so, we would infer a breakdown scale |q| = 600 MeV, 
as there the short-distance effect ~ L2 becomes equal in 
magnitude to the impulse-approximation result. 



Order 


A/ A [MeV] 


dg [GeV-2] 


L2 [GeV-"] 


NLO 


400/500 


-0.010 


0.243 


NLO 


400/700 


-0.011 


0.249 


NLO 


550/500 


0.016 


0.605 


NLO 


550/600 


0.017 


0.731 


NLO 


550/700 


0.018 


0.892 


NNLO 


450/500 


-0.011 


0.188 


NNLO 


450/700 


-0.009 


0.173 


NNLO 


550/600 


0.005 


0.089 


NNLO 


600/500 


0.001 


0.113 


NNLO 


600/700 


-0.001 


0.028 



TABLE L Values for dg and L2 found by fitting data up to 
|q| = 400 MeV, using different values of the cutoffs. 
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In Table U we present the values of dg and L2 obtained 
in our fits. Small values of dg a re p referred, which is 
consistent with the findings of Ref. [23 • Reassuringly, the 
inferred values of dg show only a very mild dependence 
on the cutoffs as compared to the expected natural size of 
this LEG, |(ig| ^ 1 GcV^^. In contrast, the values of L2 
do depend on the choice of the regulator employed for the 
NN potential, as one would expect. It is comforting to 
see that all obtained values of L2 are natural with respect 
to the cutoff scale A employed in these calculations. The 
values of L2 reported in the table show that two-body 
effects in Gm play a larger role in the calculation with 
NLO deuteron wave functions, as seen in Fig. [2] 

Summary: The first two-body effects in the deuteron 
magnetic form factor Gm, occur at 0{eP'^) in xEFT, 
i.e. three orders beyond leading. Inclusion of these 
mechanisms in the computation of Gm improves the de- 
scription of data, and allows exact reproduction of the 
deuteron magnetic moment, which otherwise is underpre- 
dicted in xEFT. Experimental data is then well described 
for < 0.35 GeV^, and the chiral expansion for Gm is 



found to converge well for < 0.25 GeV^, provided 
that the NNLO wave functions of Ref. Q are employed. 
Finally, we note that the proposal of Ref. ^ to change 
the scaling of short-distance xEFT operators in order to 
ensure proper renormalization of the theory does not sig- 
nificantly alter the relative importance of such operators 
in the ^Si-'^Di channel Therefore we expect that the 
conclusions of this study will be quite robust with respect 
to developments on this front. 
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